4.1 


INTEGRATION 


THE DEFINITE INTEGRAL 


We shall begin our study of the integral calculus in the same way in which we began 
with the differential calculus—by asking a question about curves in the plane. 
Suppose f is a real function continuous on an interval I and consider the 
curve y = f(x). Let a < b where a, b are two points in J, and let the curve be above the 
x-axis for x between a and b; that is, f(x) => 0. We then ask: What is meant by the 
area of the region bounded by the curve y = f(x), the x-axis, and the lines x = a and 
= b? That is, what is meant by the area of the shaded region in Figure 4.1.1? We 
call this region the region under the curve y = f(x) between a and b. 


Figure 4.1.1 The Region under a Curve 


The simplest possible case is where f is a constant function; that is, the curve 
is a horizontal line f(x) = k, where k is a constant and k > 0, shown in Figure 4.1.2. 
In this case the region under the curve is just a rectangle with height k and width 
b — a, so the area is defined as 


Area = k+(b — a). 


The areas of certain other simple regions, such as triangles, trapezoids, and semi- 
circles, are given by formulas from plane geometry. 
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area = k(b — a) 
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Figure 4.1.2 


The area under any continuous curve y = f(x) will be given by the definite 
integral, which is written 


fy (x) dx. 


Before plunging into the detailed definition of the integral, we outline the main ideas. 

First, the region under the curve is divided into infinitely many vertical 
strips of infinitesimal width dx. Next, each vertical strip is replaced by a vertical 
rectangle of height f(x), base dx, and area f(x) dx. The next step is to form the sum 
of the areas of all these rectangles, called the infinite Riemann sum (look ahead to 
Figures 4.1.3 and 4.1.11). Finally, the integral [° f(x) dx is defined as the standard 
part of the infinite Riemann sum. 

The infinite Riemann sum, being a sum of rectangles, has an infinitesimal 
error. This error is removed by taking the standard part to form the integral. 

It is often difficult to compute an infinite Riemann sum, since it is a sum of 
infinitely many infinitesimal rectangles. We shall first study finite Riemann sums, 
which can easily be computed on a hand calculator. 

Suppose we slice the region under the curve between a and 5 into thin vertical 
strips of equal width. If there are n slices, each slice will have width Ax = (b — a)/n. 
The interval [a, b] will be partitioned into n subintervals 


(Xo, x1], [x1 ’ Xa], aoe ee) [x3 1s Xils 
where Xp =4,X, =a+ Ax,x, =a+4+2Ax,...,x, = b. 
The points x9,X,,...,X, are called partition points. On each subinterval {x,_,.Xx,], 


we form the rectangle of height /(x,_,). The kth rectangle will have area 
P(X_- 1) + Ax. 


From Figure 4.1.3, we can see that the sum of the areas of all these rectangles will be 
fairly close to the area under the curve. This sum is called a Riemann sum and is equal 
to 


TP (Xo) Ax + f(x) Ax 4.°°°> 4+ fx, 1) Ax. 
It is the area of the shaded region in the picture. A convenient way of writing Riemann 


sums is the “X-notation” (< is the capital Greek letter sigma), 


Dd f(x) Ax = P(X) Ax + f(xy) Ax ++ + f(%,—1) Ax. 
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f() 


j 
Vs 
Figure 4.1.3. The Riemann Sum 


The a and b indicate that the first subinterval begins at a and the last subinterval ends 
at b. 

We can carry out the same process even when the subinterval length Ax does 
not divide evenly into the interval length b — a. But then, as Figure 4.1.4 shows, there 
will be a remainder left over at the end of the interval [a, b], and the Riemann sum will 
have an extra rectangle whose width is this remainder. We let n be the largest integer 
such that 


a+nAx <b, 
and we consider the subintervals 
(Xo, x1], baa J (Xn 1? Xn, {Xn bj, 
where the partition points are 


Xp =a, X,y =a+ AX, Xp =aA42Ax,..., xX, =a+ndAx, b. 


Figure 4.1.4 
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x, Will be less than or equal to b but x, + Ax will be greater than b. Then we define 
the Riemann sum to be the sum 


Sava TALE AA 2 TN Sed 


Thus given the function f, the interval [a,b], and the real number Ax > 0, we have 
defined the Riemann sum }”? f(x) Ax. We repeat the definition more concisely. 


DEFINITION 


Let a <b and let Ax be a positive real number. Then the Riemann sum 
by a . 
> f(x) Ax is defined as the sum 
b 


DS (x) Ax = f(x) Ax + S061) Ax + + fH q-1) AX + FOB = Xy) 


a 
where n is the largest integer such that a + n Ax < b, and 
Xp =a, X,=at+ Ax,'', x,=a+tnAx, b 


are the partition points. 


If x, = b, the last term f(x,)(b — x,) is zero. The Riemann sum 5” f(x) Ax 
is a real function of three variables a, b, and Ax, 
BD 
¥ f(x) Ax = S(a, b, Ax). 
The symbol x which appears in the expression is called a dummy variable (or bound 
variable), because the value of $°? f(x) Ax does not depend on x. The dummy variable 
allows us to use more compact notation, writing f(x) Ax just once instead of writing 
I (Xo) Ax, f(x) Ax, f(x) Ax, and so on. 
From Figure 4.1.5 it is plausible that by making Ax smaller we can get the 
Riemann sum as close to the area as we wish. 


F(x) 


Figure 4.1.5 


EXAMPLE 1 Let f(x) = 4x. In Figure 4.1.6, the region under the curve from x = 0 
to x = 2 is a triangle with base 2 and height 1, so its area should be 
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Figure 4.1.6 


Let us compare this value for the area with some Riemann sums. In Figure 
4.1.7, we take Ax = 4. The interval [0,2] divides into four subintervals 
(0, 4], [4 1], (1, 31, and [, 2]. We make a table of values of f (x) at the lower 


endpoints. 
Xi 0441 2 
FOS |) 02 4 


I 
Ax= 7 


Riemann sum = $ 


Figure 4.1.7 


The Riemann sum is then 


In Figure 4.1.8, we take Ax = 4. The table of values is as follows. 
my [0422344 § 3 
fx) |04%3 ¢@ # € 3 8 3 


The Riemann sum is 


2 
Vf@)Ax=O-f+ SEF RStS at Sst Fat Seat gs Hs 
0 


We see that the value is getting closer to one. 


Finally, let us take a value of Ax that does not divide evenly into the interval 
length 2. Let Ax = 0.6. We see in Figure 4.1.9 that the interval then divides 
into three subintervals of length 0.6 and one of length 0.2, namely [0, 0.6], 
[0.6, 1.2], [1.2, 1.8], (1.8, 2.0]. 


X, 0 06 12 18 
f(x, | 0 03 06 09 
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| Y| Ax = 0.6 
Riemann sum = .72 


Ax = 7 


Riemann sum = 


aX 
\\ 


KQY 


O 


Figure 4.1.8 Figure 4.1.9 
The Riemann sum is 


2 
Yd f(x) Ax = 0(.6) + (.3)(.6) + (.6)(.6) + (.9)(.2) = .72. 
0 


EXAMPLE 2 Let f(x) = \/1 — x’, defined on the closed interval J = [—1, 1]. The 
region under the curve is a semicircle of radius 1. We know from plane 
geometry that the area is 2/2, or approximately 3.14/2 = 1.57. Let us compute 
the values of some Riemann sums for this function to see how close they are 
to 1.57. First take Ax = 3 as in Figure 4.1.10(a). We make a table of values. 


x =P 12 0 172 
fe) | 0 (34 1 /3/4 


The Riemann sum is then 
1 
Y f(x) Ax = 0+1/2 + /3/4-1/2 + 1+ 1/2 + /3/4- 1/2 
=1 


= 1+/3 a 1,97. 
2 


Next we take Ax = 4. Then the interval [ —1, 1] is divided into ten subinter- 
vals as in Figure 4.1.10(b). Our table of values is as follows. 


i, Pag ot ee 2 alg ok Ese 
| a ee SP 2 BBS 
fol a ee a5 AES 

a a een She 8 


Yi 


Figure 4.1.10 
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The Riemann sum is 


. 1 3 4 21 ./24 
O+z+e¢+ + 
5 5 5 5 5 


19 + 2/21 + 2,/24 es 
= 5 ieee Ue . 


JB. Ji 


+1+ + gee ie 
5 5. aor 


— 


Thus we are getting closer to the actual area n/2 ~ 1.57. 


By taking Ax small we can get the Riemann sum to be as close to the area 
as we wish. 

Our next step is to take Ax to be infinitely small and have an infinite Riemann 
sum. How can we do this? We observe that if the real numbers a and 3b are held fixed, 
then the Riemann sum 


¥ feo Ax = S(Ax) 


is a real function of the single variable Ax. (The symbol x which appears in the 
expression is a dummy variable, and the value of 


) J (x) Ax 
depends only on Ax and not on x.) Furthermore, the term 
fe) Ax = S(Ax) 
is defined for all real Ax > 0. Therefore by the Transfer Principle, 
¥ fo) dx = S(dx) 


is defined for all hyperreal dx > 0. When dx > 0 is infinitesimal, there are infinitely 
many subintervals of length dx, and we call 


b 
Y fs) dx 


an infinite Riemann sum (Figure 4.1.11). 


Figure 4.1.11 Infinite Riemann Sum 
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We may think intuitively of the Riemann sum 


b 
¥ f(x) ax 
as the infinite sum 


Sf (Xo) dx + f(x)dx +++ + flxq_1) dx + f(xq)(b — Xx) 


where H is the greatest hyperinteger such that a + Hdx < b. (Hyperintegers are 
discussed in Section 3.8.) H is positive infinite, and there are H + 2 partition points 
Xo, Xj,+.+;Xq, b. A typical term in this sum is the infinitely small quantity f(x) dx 
where K is a hyperinteger,0 < K < H,andxx =a+ K dx. 

The infinite Riemann sum is a hyperreal number. We would next like to take 
the standard part of it. But first we must show that it is a finite hyperreal number and 
thus has a standard part. 


THEOREM 1 


Let f be a continuous function on an interval I, let a < b be two points in I, and 
let dx be a positive infinitesimal. Then the infinite Riemann sum 


b 
¥ f(x) dx 
is a finite hyperreal number. 


PROOF Let B be a real number greater than the maximum value of f on [a, b]. 
Consider first a real number Ax > 0. We can see from Figure 4.1.12 that the 


| 
7) 
ma 


Figure 4.1.12 a b 


finite Riemann sum is less than the rectangular area B«(b — a); 
eye < B+(b — a). 

Therefore by the Transfer Principle, 
x f(x) dx < B+(b — a). 


In a similar way we let C be less than the minimum of f on [a,b] and show 
that 
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b 
¥ f(x)dx > C+(b — a). 


Thus the Riemann sum $? f(x)dx is finite. 


We are now ready to define the central concept of this chapter, the definite 
integral. Recall that the derivative was defined as the standard part of the quotient 
Ay/Ax and was written dy/dx. The “definite integral” will be defined as the standard 
part of the infinite Riemann sum 

b 


> f(x) dx, 
and is written j° f(x) dx. Thus the Ax is changed to dx in analogy with our differential 
notation. The & is changed to the long thin S, i.e., f, to remind us that the integral is 
obtained from an infinite sum. We now state the definition carefully. 


DEFINITION 


Let f be a continuous function on an interval I and let a < b be two points in I. 
Let dx be a positive infinitesimal. Then the definite integral of f from a to b with 
respect to dx is defined to be the standard part of the infinite Riemann sum with 
respect to dx, in symbols 


[109 dx = af SF) dx 


We also define | T(x) dx = 0, 


i oravde= - [709 ae 


By this definition, for each positive infinitesimal dx the definite integral 


ie: (x) dx 


is a real function of two variables defined for all pairs (u, w) of elements of J. The 
symbol x is a dummy variable since the value of 


[ f(x) dx 


does not depend on x. 

In the notation }”? f(x) dx for the Riemann sum and f° f(x) dx for the 
integral, we always use matching symbols for the infinitesimal dx and the dummy 
variable x. Thus when there are two or more variables we can tell which one is the 
dummy variable in an integral. For example, x7t can be integrated from 0 to 1 with 
respect to either x or ¢. With respect to x, 


1 
dx? dx = xptdx + xftdx +++: + x#_,tdx 
0 
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(where dx = 1/H), and we shall see later that 
1 
| x?tdx = st(xétdx + x?tdx +++: + x},_,tdx) = 40. 
0 


With respect to t, however, 


I 
Sattar = xt tg dt + x7h dt hae x7 ty sdt 
0 


and we shall see later that 
1 
i x?*rdt = 4x?. 
0 


The next two examples evaluate the simplest definite integrals. These 
examples do it the hard way. A much better method will be developed in Section 4.2. 


EXAMPLE 3 Given a constant c > 0, evaluate the integral |? c dx. 


Figure 4.1.13 shows that for every positive real number Ax, the finite Riemann 
sum is 


b 
Yc Ax = c(b — a). 


By the Transfer Principle, the infinite Riemann sum in Figure 4.1.14 has the 
same value, 


b 
Yi cdx = c(b — a). 
Taking standard parts, 


b 
[eax = c(b — a). 


a 


This is the familiar formula for the area of a rectangle. 


Figure 4.1.13 Figure 4.1.14 
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EXAMPLE 4 Given b > 0, evaluate the integral § x dx. 


The area under the line y = x is divided into vertical strips of width dx. 
Study Figure 4.1.15. The area of the lower region A is the infinite Riemann 
sum 


b 
(1) area of A= ) x dx. 
0 


By symmetry, the upper region B has the same area as A; 
(2) area of A = area of B. 


Call the remaining region C, formed by the infinitesimal squares along the 
diagonal. Thus 


(3) area of A + area of B + area of C = b?. 


Each square in C has height dx except the last one, which may be smaller, 
and the widths add up to b, so 


(4) 0 < area of C < bdx. 
Putting (1)-(4) together, 


b b 
2). kde sb? < (2 xas} + bdx. 
0 0 


Since b dx is infinitesimal, 
5 
23) x dx = b, 
0 


b b2 


Taking standard parts, we have 


Figure 4.1.15 
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PROBLEMS FOR SECTION 4.1 


Compute the following finite Riemann sums. If a hand calculator is available, the Riemann sums 
can also be computed with Ax = 75. 


1 yi (3x + 1) Ax, Ax =4 2 Yo Bx + I Ax, Ax =3 
3 Yt, Bx + 1)Ax, Ax=$ 4 Yip 2x? Ax, Ax = 
5 Si yr Ax, Ax =1 6 ¥6 (2x — 1) Ax, Ax =1 
a Ys (2x — 1) Ax, Ax =2 8 Ye 1 ? — Ax, Ax = 
9 Foot Sax, Ax=4 10 gle Sas, Ax = ¥5 
a ¥3 (5x? — 12)Ax, Ax = 2 12 ¥3 (5x? — 12)Ax, Ax =1 
13 Did + 1/x) Ax, Ax =4 14 VF 10" Ax: Ax =4 
15 yo ot ax, Ax =4 16 yy 2k As Ax =} 
17 pay ace Ax = 1 18 D2 x — 4 Ax, Ax = 2 
19 % sinx Ax, Ax=n4 20 Ya sin? x Ax. Ax = 714 
21 yh et Ax, Ax=1/5 22 yh xe* Ax, Ax = 1/5 
23 Yi} nx Ax, Ax = 1 24 ye mn Ax. Ax =1 
25 Let b be a positive real number and n a positive integer. Prove that if Ax = b/n, 
b 
yx Ax = (142 4°+++ (2 — 1) Ax?. 
0 
ee a(n — 1) 
Using the formula 1 + 24+ °::+(n-— I= 7? Prove that 
b 
yx Ax = (1 = L/n)b?/2. 
O 
26 Let H bea positive infinite hyperinteger and dx = b/H. Using the Transfer Principle and 
Problem 25, prove that J} x dx = b?/2. 
27 Let b be a positive real number, n a positive integer, and Ax = b/n. Using the formula 
—1 -1 
127427437 4°54 (n- y= oe, 
prove that 
u — 1)(2n — 1) b3 
Bik n(n 
2» ‘ 6 ” 
28 Use Problem 27 to show that [3 x* dx = b°/3. 


FUNDAMENTAL THEOREM OF CALCULUS 


In this section we shall state five basic theorems about the integral, culminating in 
the Fundamental Theorem of Calculus. Right now we can only approximate a 
definite integral by the laborious computation of a finite Riemann sum. At the end 
of this section we will be in a position easily to compute exact values for many definite 
integrals. The key to the method is the Fundamental Theorem. Our first theorem 
shows that we are free to choose any positive infinitesimal we wish for dx in the 
definite integral. 


4.2 FUNDAMENTAL THEOREM OF CALCULUS 


THEOREM 1 


Given a continuous function f on [a,b] and two positive infinitesimals dx 
and du, the definite integrals with respect to dx and du are the same, 


[109 dx = [re du. 


From now on when we write a definite integral [? f(x) dx, it is understood 
that dx is a positive infinitesimal. By Theorem 1, it doesn’t matter which infinitesimal. 

The proof of Theorem 1 is based on the following intuitive idea. Figure 4.2.1 
shows the two Riemann sums )? f(x) dx and )? f(u) du. We see from the figure 
that the difference )® f(x) dx — )® f(u) du is a sum of rectangles of infinitesimal 
height. These difference rectangles all lie between the horizontal lines y = —é and 
y = e, where ¢ is the largest height. Thus —e(b — a) < )? f(x)dx — ¥2 f(wdu < 
e(b — a). Taking standard parts, 


b 
O< Preoax - i) f@du <0, 


[ oa { ” pu) du. 


Figure 4.2.1 
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Theorem 1 shows that whenever Ax is positive infinitesimal, the Riemann 
sum is infinitely close to the definite integral, 


b b 
ye Sf (x) Ax | f(x) dx. 


This fact can also be expressed in terms of limits. It shows that the Riemann sum 
approaches the definite integral as Ax approaches 0 from above, in symbols 


b 
ff f(x) dx = lim f(x) Ax. 


Ax7>Ot a4 


Given a continuous function f on an interval I, Theorem 1 shows that the 
definite integral is a real function of two variables a and b, 


b 
A(a, b) = S(x) dx, a, bin I, 


We now formally define the area as the definite integral shown in Figure 4.2.2, 


Figure 4.2.2 


DEFINITION 


[If f is continuous and f(x) = 0 on [a,b], the area of the region below the 
curve y = f(x) from a to b is defined as the definite integral: 


b 
Area = i Sx) dx. 


The next two theorems give basic properties of the integral. 


THEOREM 2 (The Rectangle Property) 


Suppose f is continuous and has minimum value m and maximum value M 
on a Closed interval [a,b]. Then 


b 
mb—-—a)< { S(x)dx < Mb — a). 


That is, the area of the region under the curve is between the area of the rectangle 
whose height is the minimum value of f and the area of the rectangle whose 
height is the maximum value of f in the interval [a, b]. 


4.2 FUNDAMENTAL THEOREM OF CALCULUS 


The Extreme Value Theorem is needed to show that the minimum value m 
and maximum value M exist. The rectangle of height m is called the inscribed rectangle 
of the region, and the rectangle of height M is called the circumscribed rectangle. 
From Figure 4.2.3, we see that the inscribed rectangle is a subset of the region under 
the curve, which is in turn a subset of the circumscribed rectangle. The Rectangle 
Property says that the area of the region is between the areas of the inscribed and 


MMI. 
b 


a 


Figure 4.2.3 The Rectangle Property 
PROOF By Theorem 1, any positive infinitesimal may be chosen for dx. Let us 
choose a positive infinite hyperinteger H and let dx = (b — a)/H. Then 


dx evenly divides b — a; that is, the interval [a,b] is divided into H sub- 
intervals of exactly the same length dx. Then 


b 
Y mdx = m+H+dx = m(b — a), 
b 
YM dx = M+ H+dx = M(b — a). 


For each x, we have m < f(x) < M. Adding up and taking standard parts, 
we obtain the required formula. 


b b b 
y mdx <¥ f(x)dx < ¥ Max, 


m(b—a)< [ fooax < M(b — a). 


One useful consequence of the Rectangle Property is that the integral of 
a positive function is positive and the integral of a negative function is negative: 


b 

If f(x) > 0 on [a, 5], then 0 < m(b — a) < f(x) dx. 
b 

If f(x) < 0 on [a, 5], then [ f(x)dx < M(b — a) <0. 


The definite integral of a negative function f(x) = —g(x) from a to } is 
just the negative of the area of the region above the curve and below the x axis. 
This is because 


f(x) dx = —g(x) dx, 
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b b 
Y fi) dx = — ¥ g(x) dx, 
_ a 

} f()dx = — | g(x) dx. 


(See Figure 4.2.4.) 


Figure 4.2.4 


THEOREM 3 (The Addition Property) 


Suppose f is continuous on an interval I. Then for all a, b,c in I, 


[ 109 ax= [ 409 dx + [100 dx. 


This property is illustrated in Figure 4.2.5 for the case a<b<c. The 
Addition Property holds even if the points a,b,c are in some other order on the 
real line, such asc <a < b. 


Figure 4.2.5 
PROOF First suppose that a < b < c. Choose a dx that evenly divides the first 


interval length b — a. This simplifies our computation because it makes 
ba partition point, b = a + H dx. Then, as Figure 4.2.6 suggests, 


c b ¢ 
Y f@xydx = ¥ f(x)dx + ¥ f(x) dx. 
a a b 

Taking standard parts we have the desired formula 


[10 dx = [ 409 dx + [10 dx. 


4.2 FUNDAMENTAL THEOREM OF CALCULUS 


LZ] 


Yip 


Z| 
AY 
OO 


Figure 4.2.6 a c 


To illustrate the other cases, we prove the Addition Property when 
c <a <b. The previous case gives 


[seve = [fooax + [reve 


Since reversing the endpoints changes the sign of the integral, 


c c b 
_ dx = — d dx, 
[ 4 x i f(x) dx + I f(xax 
and the desired formula 
c b c 
[ sepa = i S(x) dx + [ reoax 
a a b 
follows. 


The definite integral of a curve can be thought of as area even if the curve 
crosses the x-axis. The curve in Figure 4.2.7 is positive from a to b and negative from 
b to c, crossing the x-axis at b. The integral {° f(x) dx is a positive number and the 
integral {§ f(x)dx is a negative number. By the Addition Property, the integral 


[ fooas = [ seax + [ seoas 


is equal to the area from a to b minus the area from b to c. The definite integral 
J. f(x) dx always gives the net area between the x-axis and the curve, counting 
areas above the x-axis as positive and areas below the x-axis as negative. 

The definite integral [% f (2) dt is a real function of two variables u and v 
and does not depend on the dummy variable t. If we replace u by a constant a and v 
by the variable x, we obtain a real function of one variable x, given by 


F(x) = i" f(t) dt. 


Our fourth theorem states that this new function is continuous. 


Figure 4.2.7 
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THEOREM 4 


Let f be continuous on an interval I. Choose a point a in I. Then the function 
F(x) defined by: 


F(x) = i “Sf (tat 
is continuous on I. 


SKETCH OF PROOF Let c be in J, and let x be infinitely close to c and between 
the endpoints of /. By the Addition Property, 


i f(t = [ “fd + i: fat, 
[10 dt — [ 10 dt = [so dt, 


and F(c) — F(x) = [ 10 dt. 


This is the area of the infinitely thin strip under the curve » = f(t) between 
t= x and t = c (see Figure 4.2.8). The strip has width Ax = c — x. By the 
Rectangle Property, its area is between m Ax and M Ax and hence is infinitely 
small. Therefore F(x) is infinitely close to F(c), and F is continuous on J. 


AW 


FO ~ FO 


b 


Figure 4.2.8 a c 


Our fifth theorem, the Fundamental Theorem of Calculus, shows -that 
the definite integral can be evaluated by means of antiderivatives. The process of 
antidifferentiation is just the opposite of differentiation. To keep things simple, let J 
be an open interval, and assume that all functions mentioned have domain /. 


DEFINITION 


Let f and F be functions with domain I. If f is the derivative of F, then F is 
called an antiderivative of f. 


4.2 FUNDAMENTAL THEOREM OF CALCULUS 


For example, suppose a particle is moving upward along the y-axis with 
velocity v = f(t) and position y = F(t) at time t. The position y = F(t) is an anti- 
derivative of the velocity v = f(t). We shall discuss antiderivatives in more detail 
in the next section. We are now ready for the Fundamental Theorem. 


FUNDAMENTAL THEOREM OF CALCULUS 


Suppose f is continuous on its domain, which is an open interval I. 


(i) For each point a in I, the definite integral of f from a to x considered as a 
function of x is an antiderivative of f. That is, 


af soa) = f(x) dx. 


(ii) If F is any antiderivative of f, then for any two points (a,b) in I the 
definite integral of f from a to b is equal to the difference F(b) — F(a), 


[169 dx = F(b) — F(a). 


The Fundamental Theorem of Calculus is important for two reasons. First, 
it shows the relation between the two main notions of calculus: the derivative, which 
corresponds to velocity, and the integral, which corresponds to area. It shows that 
differentiation and integration are “inverse” processes. Second, it gives a simple 
method for computing many definite integrals. 


EXAMPLE 1 

(a) Find fe c dx. Since cx is an antiderivative of c, 

b 

| cdx = cb ~ ca = c(b — a). 

a 
(b) Find f” x dx. 3x? is an antiderivative of x. Thus 

6 
i) x dx = 5b? — 4a’. 


The above example gives the same result that we got before but is much 
simpler. We can easily go further. 


EXAMPLE 2 Find J? x? dx. x*/3 is an antiderivative of x? because 


d(x3/3) _ 3x? 3 


dx ee 
b b3 3 
Therefore [ x? dx = a a 


This gives the area of the region under the curve y = x? between a and b 
(Figure 4.2.9). 
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Figure 4.2.9 


If a particle moves along the y-axis with continuous velocity v = f(r), the 
position y = F(t) is an antiderivative of the velocity, because v = dy/dt. The 
Fundamental Theorem of Calculus shows that the distance moved (the change in }) 
between times t = a and t = b is equal to the definite integral of the velocity, 


b 
distance moved = F(b) — F(a) = i fo dt. 


EXAMPLE 3 A particle moves along the y-axis with velocity v = 8t* cm/sec. How 


far does it move between times t= —I1 and t = 2sec? The function 
G(t) = 2r* is an antiderivative of the velocity v = 8t*. Thus the definite 
integral is 


2 
distance moved = | 823 dt = 2°24 —2+(—1)* = 30cm. 
al 


EXAMPLE 4 Find {§ Jt dt (Figure 4.2.10). The function Jt. is defined and con- 
tinuous on the half-open interval [0, oc). But to apply the Fundamental 
Theorem we need a function continuous on an open interval that contains 
the limit points 0 and 4. We therefore define 


0 fort <0 
9 = Ji for t > 0. 


This function is continuous on the whole real line. In particular it is con- 
tinuous at 0 because if t ~ 0 then f(t) = 0. The function 


Fit) = " fort <0 


232 fort > 0 


Figure 4.2.10 
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is an antiderivative of f. Then 


4 
| Jtdt = F(4) — F(0) = §+ 49? — 3.07) = 38, 
0 


In the next section we shall develop some methods for finding antiderivatives, 
The antiderivative of a very simple function may turn out to be a “new” function 
which we have not yet given a name. 


EXAMPLE 5 The only way we can show that the function f(x) = ./1 + x* has an 
antiderivative is to take a definite integral 


ee + t* dt. 
0 


This is a “new” function that cannot be expressed in terms of algebraic, 
trigonometric, and exponential functions without calculus. 


The Fundamental Theorem can also be used to find the derivative of a 
function which is defined as a definite integral with a variable limit of integration. 
This can be done without actually evaluating the integral. 


2 x 
EXAMPLE 6 Let y = [vi + t? dt. Then y = -| J1+ tdi, 
x 2 
and ay = ~af VIFF a = - JT FH ax 
2 


1 
= —— dF. 
Past 
Let u = x? + x, Then 


wy d 1 
GD Gyo if, y={ ey 
dx 3 


By the Chain Rule, 


dy dydu 1 2x +1 
= = A Dae as, 
dx dudx malo. (x? +x) 41 


x?+x 
EXAMPLE 7 Let y -{ 
3 


We conclude this section with a proof of the Fundamental Theorem of 
Calculus. 


PROOF (i) Let F(x) be the area under the curve y = f(t) from a to x, 
F(x) = | f@ dt. 


Imagine that the vertical line cutting the t-axis at x moves to the right as 
in Figure 4.2.11. 
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Figure 4.2.11 


We show that the rate of change of F(x) is equal to the length f(x) of the 
moving vertical line. 


Suppose x increases by an infinitesimal amount Ax > 0. Then 
xtAx 
F(x + Ax) — F(x) = { f(t)dt 


is the area of an infinitely thin strip of width Ax and height infinitely close to 
f(x). By the Rectangle Property the area of the strip is between the inscribed 
and circumscribed rectangles (Figure 4.2.12), 


mAx < F(x + Ax) — F(x) < M Ax. 
F(x + Ax) — F(x) 


Dividing by Ax, m < ———_———— < M. 
Ax 


Since f is continuous at x, the values m and M are both infinitely close to f(x), 
and therefore 
F(x + Ax) — F 
Ax 


Oy f(x). 


The proof is similar when Ax < 0. Hence F’(x) = f(x). 


F(xt+ Ax) - F(x) 


Figure 4.2.12 


4.2 FUNDAMENTAL THEOREM OF CALCULUS 
PROOF (ii) Let F(x) be any antiderivative of f/ Then, by (), 
a{ Foo — f sear) = Fo) - Fos) =0. 


In Section 3.7 on curve sketching, we saw that every function with derivative 
zero is constant. Thus 


F(x) - [ f@dt =Co, F(x) = [ro dt+ Cp 
for some constant Cy. Then 
b a 
F(b) — F(a) = ({ f@dt + Ca] — ({ S(t)}dt + C6) 
b b 
= [ s@ae—o= | soar 


b 
so F(b) — F(a) = | f(x) dx. 


PROBLEMS FOR SECTION 4.2 


In Problems 1-14, find an antiderivative of the given function. 


1 f(x) = 8\/x 2 f(x) = 4//x 
3 f®=3? +1 4 f(x) = 5x3 
5 f@) = 4 —- 30 6 f(z) = 2/2? 
t f(s) = 7s~° 8 f= 40? 
9 f(x) = (x — 6? 10 f(u) = (Su + 1? 
1 fy) = y? 12 f(x) = 2/x,/x 
13 fd) = |xl 14 F(t) = [2t - 4| 
15 If F(x) = x + x? for all x, find F(1) — F(—1). 
16 If F(x) = x* for all x, find F(2) — F(1). 
17 If F(é) = t*? for all t, find F(8) — F(0). 
Evaluate the definite integrals in Problems 18-22. 
1 2 
18 2x? dx 19 i x3 dx 
-t -2 
Si 4. 
20 { t dt 21 i 2./x dx 
-2 18) 
-2 
22 i —5x* dx 
-3 


In Problems 23-27 an object moves along the y-axis. Given the velocity v, find how far the object 
moves between the given times fy and f,. 


23 v=2t4+ 5, tj =0, t, =2 
24 v=4-t, f=1, t,=4 
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25 
26 
27 


28 
29 
30 
31 
32 
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ree ia 2. 1 S6 
t= 3¢?, St are 3 
v= 10077, fo= 1, t, = 100 


y= 4— x’, a=-2, b=2 
pen / e+ 2; a=-2, b=2 
y = 9x — x, a=0, b=3 
Peery a=0, b=1 
y = 3x3, a=1, b=8 


If F(t) = +t — 1 for all rand F(0) = 2, find F(2). 

If F(x) = 1 — x? for all x and F(3) = 5, find F(—1). 

Suppose F(x) and G(x) have continuous derivatives and F(x) + G(x) = 0 for all x. 
Prove that F(x) + G(x) is constant. 


Suppose F(x) and G(x) have continuous derivatives such that F(x) < G(x) for all x. 
Prove that F(b) — F(a) < G(b) — Gla) 

where a < b. 

Prove that a function F(x) has a constant derivative if and only if F(x) is linear, i-e., of the 
form F(x) = ax + b. 

Prove that a function F(x) has a constant second derivative if and only if F(x) has the 
form F(x) = ax?+ bx +. 

Suppose that F”(x) = G”(x) for all x. Prove that F(x) and G(x) differ by a linear function, 
that is, G(x) = F(x) + ax + b for some real numbers a and b. 


4.3 INDEFINITE INTEGRALS 


The Fundamental Theorem of Calculus shows that every continuous function f 
has at least one antiderivative, namely F(x) = J? f(t) dt. Actually, f has infinitely 
many antiderivatives, but any two antiderivatives of f differ only by a constant. This 
is an important fact about antiderivatives, which we state as a theorem. 


THEOREM 1 


Let f be a real function whose domain is an open interval I. 


QQ) If F(x) is an antiderivative of f(x), then F(x) + C is an antiderivative 
of { (x) for every real number C. 

(i) If F(x) and G(x) are two antiderivatives of f(x), then F(x) — G(x) is 
constant for all x in I. That is, 


G(x) = F(x) + € 


for some real number C. 


Discussion Parts (i) and (ii) together show that if we can find one antiderivative 


F(x) of f(x), then the family of functions 
F(x) + C, C = areal number 


4.3 INDEFINITE INTEGRALS 199 


gives all antiderivatives of f(x). We see from Figure 4.3.1 that the graph 
of F(x) + C is just the graph of F(x) moved vertically by a distance C. The 
graphs of F(x) and F(x) + C have the same slopes at every point x. For 
example, let f(x) = 3x?. Then F(x) = x? is an antiderivative of 3x? because 


d(x*) | 
dx 
But x° + 6 and x? — af are also antiderivatives of 3x7. In fact, x* + C is 


an antiderivative of 3x? for each real number C. Theorem 1 shows that 3x? 
has no other antiderivatives. 


3x?, 


Figure 4.3.1 


PROOF We prove (i) by differentiating, 


d(F(x) + C) x d(F(x)) | dc 


i. Hy + Fe =I) + 0 =f). 


Part (ii) follows from a theorem in Section 3.7 on curve sketching. If a 
function has derivative zero on I, then the function is constant on J. The 
difference F(x) — G(x) has derivative f(x) —f(x) =0 and is therefore 
constant. We used this fact in the proof of the Fundamental Theorem of 
Calculus. 


In computing integrals of f, we usually work with the family of all anti- 
derivatives of f. We shall call this whole family of functions the indefinite integral of f- 
The symbol for the indefinite integral is | f(x) dx. If F(x) is one antiderivative of f, 


the indefinite integral is the set of all functions of the form F(x) + Cy, Cp constant. 
We express this with the equation 


[ feo ax =Fx)+C. 


It is an equation between two families of functions rather than between two single 
functions. C is called the constant of integration. To illustrate the notation, 


[ 3x2 ax = x° AG. 


We repeat the above definitions in concise form. 
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DEFINITION 


Let the domain of f be an open interval I and suppose f has an antiderivative. 
The family of all antiderivatives of f is called the indefinite integral of f and is 
denoted by { f(x) dx. 


Given a function F, the family of all functions which differ from F only by a 
constant is written F(x) + C. Thus if F is an antiderivative of f we write 


[ fooax = F(x) + C. 


When working with indefinite integrals, it is convenient to use differentials 
and dependent variables. If we introduce the dependent variable u by u = F(x), then 


du = F(x) dx = f(x) dx. 


Thus the equation [709 dx = F(x) + C 


can be written in the form [aw =u+C. 


The differential symbol d and the indefinite integral symbol | behave as 
inverses to each other. We can start with the family of functions u + C, form du, and 
then form | du = u + C to get back where we started. Some of the rules for differentia- 
tion given in Chapter 2 can be turned around to give a set of rules for indefinite 
integration, 


THEOREM 2 


Let u and v be functions of x whose domains are an open interval I and suppose 
du and dv exist for every x in I. 


— 


1) au =u+C. 


(ii) | Constant Rule [eau = c| du. 


(iii) Sum Rule [au +dv= [au + [a 


wer) 
(iv) Power Rule [« du = 4 + C, where r is rational, r 4 —1, 


andu>OQOonl. 


(v) | sin udu = —cosut+C, 
(vi) { c0s udu = sinu+C. 


(vii) { edu =e" 4+C. 
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1 
(viii) iF du = In|u|+C (u # 0). 


Discussion The Power Rule gives the integral of uv” whenr 4 — 1, while Rule 
(viii) gives the integral of wu” when r = —1. When we put u = f(x) and 
v = g(x), the Constant and Sum Rules take the form 


Constant Rule [ re dx = ef St (x) dx. 


Sum Rule Ker + g(x)) dx = { #09 dx + [ 00 dx. 


In the Constant and Sum Rules we are multiplying a family of functions 
by a constant and adding two families of functions. If we do either of these 
two things to families of functions differing only by a constant, we get another 
family of functions differing only by a constant. For example, 


1(3x* + C) = 21xt + 7C = 21xt + C 
is the family of all functions equal to 21x* plus a constant. Similarly, 
(3./x + C) + (Sx — «/x + D) = 5x 42/x + (C+D) =5x42/x4C 
is the family of all functions equal to 5x + 2./x plus a constant. 


PROOF OF THEOREM 2 


(i) This is just a short form of the theorem that u + C is the family of all 
functions which have the same derivative as u. 


(ii) We have c du = d(cu), whence 
[eau =f atew = cu + C= elu + Cc) = c | au 
(iii) du + dv = d(u + v), 
favsdo= faw+yautv+c=fdws fae 


r+ ” r 
d e = Ga DE du = u’ du, 
r+1 r+1 


(iv) 


[uqu-*S 4c 
r4+il 


Rules (v)—(viii) are similar. Only the last formula, (viii), requires an explana- 
tion. The absolute value in In|u| comes about by combining the two cases u > 0 
and u < 0. When u > 0, uv = [u| and 
1 
dn |u|) = din u) = aft 


When u < 0, Inu is undefined, but |u| = —u and In |u| = In (—u). Thus 


din |u|) = dan (—u)) = — * a(—w) — «du 
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Thus, in both cases, when u # 0, 


1 
d(ln |u|) = —du, 
u 
1 
[au =Injul + C. 
u 
EXAMPLE 1 [ox + 3sinx)dx = 2In|{x| —3cosx +C, 


We can use the rules to write down at once the indefinite integral of any 
polynomial. 


EXAMPLE 2 Jee — 6x? 4 2x + l)dx = x* — 2x9 4x? + x4C. 


EXAMPLE 3 ie = V5 Aes = re zx -é 


Indefinite integration is much harder than differentiation, because there are 
no rules for integrating the product or quotient of two functions. It often requires 
guesswork. The short list of rules in Theorem | will help, and as this course proceeds 
we shall add many more techniques for finding indefinite integrals. 


1+ x)?fad - l-—x 


Our rules give no hint on finding this integral. However, once the answer 
is given to us we can easily prove that it is correct by differentiating, 


EXAMPLE 4 Show that [age [HE te 
i xe 


l+x 
L—x  d{(l + x)" — x)7"?) 
dx - dx 
= ay (UC) = x) (ey a x) 
= (1+ x) — x) PPB + x) + 201 — x)] 
1 
(1+ x) = x8?" 


Here is a warning that may prevent some common mistakes. 
Warning: The integral of the product of two functions is not equal to the 
product of the integrals. The same goes for quotients. That is, 


hon fewer = [fal fas 
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fxas|[ foc 1) dx ote + s +C 


For example, 


Wrong: [ x + 1)dx = 


4 x3 
+ > + C. 
x3 x? 
Correct: al x(x + I) dx = ke +xdx= 4+ F406 
[vax 
u 
Wrong: -—dx = ; 
‘ [ vax 
For example, 
[ + 1)dx 2 
Wrong : {- a Big = = a a +C 
Vx [ Vaax 3x 
/x 3 
aes + + C. 
4 2./X 


1 
Correct: [- = dx = vs + 7 dx = 3x3/2 4 2./x +C. 


The indefinite integral can be used to solve problems of the following type. 
Given that a particle moves along the y-axis with velocity v = f(t), and that at a 
certain time t = tg its position is y = y. Find the position y as a function of t. 


EXAMPLE 5 A particle moves with velocity v = 1/t?, t > 0. At time t = 2 it is at 
position y = 1. Find the position y as a function of t. We compute 


1 1 
fea=[aa=-F+¢ 


Since dy/dt = v, y is one of the functions in the family —1/t + C. We can 
find the constant C by setting t = 2 and y = 1, 


1 
y = + G; 1 —— 
t 
Then the answer is 


! 
= --4 1h. 
ieee 


The next theorem shows that in such a problem we can always find the answer 
if we are given the position of the particle at just one point of time. 
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THEOREM 3 


Suppose the domain of f is an open interval I and f has an antiderivative. Let 
P(X93)'o) be any point with xq in I. Then f has exactly one antiderivative 
whose graph passes through P. 


PROOF Let F be any antiderivative of f Then F(x) + C is the family of all anti- 
derivatives. We show that there is exactly one value of C such that the 
function F(x) + C passes through P(x9, Yo) (Figure 4.3.2). We note that all 
of the following statements are equivalent: 


1) F(x) + C passes through P(x, }'o). 
(2) F(Xo) + C = Jo. 
3) 


P(Xo, Yo) 


Figure 4.3.2 


The Fundamental Theorem of Calculus, part (ii), may be expressed briefly 
as follows, where fis continuous on J. 


Lf J f(x) dx = F(x) + C, then 
[ feos = F(b) — F(a). 
For evaluating definite integrals we introduce the convenient notation 
ria} = F(b) — Fla). 


a 


It is read “F(x) evaluated from a to b.” 
The Constant and Sum Rules hold for definite as well as indefinite integrals: 


b b 
Constant Rule i) cf (x) dx = e| T(x) dx. 


b b b 
Sum Rule i) (f(x) + g(x) dx = | S(x)dx + [ se dx. 


The Constant Rule is shown by the computation 
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b b 
[ of (x) dx = cF(b) — cF(a) = c(F(b) — F(a)) = { T(x) dx. 
The Sum Rule is similar. 


EXAMPLE 6 Evaluate the definite integral of y = (1 + D/t? from t = 1 to r= 2 
(see Figure 4.3.3). 


2 2 
d= [ee yade 
t- ob i 
2 4 2 . t~? 2 t—i 2 
= |r a+ [or dt=-_] +—_— 
j i =) = 1h 


~ oe oe cae ae: 
~ \(=2)-4 cra) + (3- aa) -3 +37 


Thus the area under the curve y = (1 + H/t* from t = 1 tot = 2is%. 


Figure 4.3.3 


EXAMPLE 7 Find the area of the region under one arch of the curve y = sinx 
(see Figure 4.3.4). 


One arch of the sine curve is between x = 0 and x = xz. The area is the 
definite integral 


n Tt 
| sin x dx = cos +| 


0 10) 
= —cos z — (—cos0) = —(—1) — (-1) =2. 


The area is exactly 2. 


y 


y= sinx 


Figure 4.3.4 
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EXAMPLE 8 Find the area under the curve y = —2x~' from x = —Stox= —1. 
(See Figure 4.3.5.) 


The area is given by the definite integral 


-1 
| —2x7! dx. 
=§ 


First compute the indefinite integral 


i —2x7!dx = =2 [tax = —2In|x}+C. 


Now compute the definite integral. 


-1 
| —2x-'!dx = ~2 in xt | 
-5 


—2(In | —1] — In| —5]) = —2(n 1 — In 5) 
= nS 3019. 


weal. 


=5 


Figure 4.3.5 


This example illustrates the need for the absolute value in the integration rule 
[xctax = In|x|+C. 
The natural logarithm In x is undefined at x = —5 and x = —1, but In |x| is defined 


for all x # 0. The absolute value sign is put in when integrating x~' and removed 
when differentiating In |x|. 


EXAMPLE 9 In computing definite integrals one must first make sure that the 
function to be integrated is continuous on the interval. For instance, 


1 1 x7! 1 
Incorrect: i 5 dx | ak ren oe) area 

es —If_y 
This is clearly wrong because 1/x* > 0 so the area under the curve cannot be 
negative. The mistake is that 1/x? is undefined at x = 0 and hence the 
function is discontinuous at x = 0. Therefore the area under the curve and 
the definite integral 
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are undefined (Figure 4.3.6). 


F(x) 
(-1, 1) 
qi, -1) 
F(x) = x 
Figure 4.3.6 
PROBLEMS FOR SECTION 4.3 
Evaluate the following integrals. 
1 fa + 2x + 3x7) dx 2 fee ~ 6x + 9) dx 
3 fa2e — 315 + 20 + Idi 4 ic + yy"? — 4y- ay 
5 fe 410) dt 6 fees — 3y?3) dy 
7 [ex — 3) dx 8 fo — 2)(2x + 1) dx 
9 i @+ ta? de 10 { (a L/ede 
i { 5 cos x dx 12 { (sin x + cos x) dx 
1B | = tax 14 J ra ees dx 
15 fa + x7 dx 16 [ 3etax 
17 { B+ /n4-2/dat 18 | - a ds 
19 [tt wey 20 fo — x2)(1 + 4x?) dx 
21 [ee + bx + e)dx 22 J tase" + a,x* + a,x + ag) dx 
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23 


25 


27 


29 


31 
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2 
i (2x — 4x7 + x5) dx 


1 
| (1 + x? + 3x*) dx 
-1 


n 
i cos x dx 


0) 
2 


1 
24 I (1 + x? + 3x*) dx 
ie) 
1 
26 { e* dx 
1. 
m2 
28 cos x dx 
0 
5 
x—1 
30 [> dx 
5) OX 


In Problems 32-36, find the position y as a function of t given the velocity vu = dy/dt and the value 
of y at one point of time. 


32 
33 
34 
35 
36 


v= 2t + 3, y=0 
v= 417 — 1, y= 2 
0 


v= 34, y= 
u = 2sint, y= 10 
pe 3k, ye] 


t=0 


when t=0 


when 


when t= —1 
when t=0 
when t=1 


In Problems 37-42, find the position y and velocity v as a function of t given the acceleration a and 
the values of y and v att = Oorr = 1. 


37 
38 
39 
40 
4] 
42 
43 


44 
45 


46 
O47 


a=t, v=0 and y=1 when tr=0 

a= —32, v=10 and y=0 when t=0 

a = 31, vc=1 and y=2 when r=0 

a=1- Jt v= -2 and y=1 when tr=0 

a=t3, e=1 and y=0 when t=1 

a= -—sint, vu=0 and j=4 when t=0 

Which of the following definite integrals are undefined? 

(a) ii mee (b) i 
ik ix 


3 
ep) | jx — I] dx 
a3 


() [ tan x dx 
0 


Find the function f such that /” is constant, {(0) = (0) and /(2) = f'(2). 


An object moves with acceleration a = 6. Find its position » as a function of f, given 
that y = | when tf = Oand y = 4 whent = I. 


Find the function ft such that h” is constant, A(1) = 1, A(2) = 2. and A(3) = 3. 


Suppose that F’(x) exists for all x. and let (xg. 49) and (x,.}';) be two given points. 
Prove that there is exactly one function G(x) such that 


4.4 INTEGRATION BY CHANGE OF VARIABLES 


G(Xo) = Yo 
G(x) =) 
G"(x) = F"(x) for all x. 
O 48 Assume that F(x) exists for all x, and let (x,, y,) and (x,, y2) be two points with x, # x). 


Prove that there is exactly one function G(x) such that G"(x) = F(x) for all x, and the 
graph of G passes through the two points (x,, y,) and (x2, y3). 


4.4 INTEGRATION BY CHANGE OF VARIABLES 


We have seen that the sum, constant, and power rules for differentiation can be turned 
around to give the sum, constant, and power rules for integration. In this section we 
shall show how to make use of the Chain Rule for differentiation in problems of 
integration. The Chain Rule will lead to the important method of integration by 
change of variables. The basic idea is to try to simplify the function to be integrated 
by changing from one independent variable to another. 

If F is an antiderivative of f and we take u as the independent variable, then 
ff) du is a family of functions of u, 


{0 du = Flu) +C. 


But if we take x as the independent variable and introduce u as a dependent variable 
u = g(x), then du and { f(u) du mean the following: 


du = g(x) dx, i f(u) du = { f(g(x)e'(x) dx = H(x) + C. 


The notation { f(u) du always stands for a family of functions of the independent 
variable, which in some cases is another variable such as x. The next theorem can be 
used as follows. To integrate a given function of x, properly choose a new variable 
u = g(x) and integrate a new function with respect to u. 


DEFINITION 


Let I and J be intervals. We say that a function g maps J into I if for every 
point x in J, g(x) is defined and belongs to I (Figure 4.4.1). 


Figure 4.4.1 g maps J into J 
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THEOREM 1 (Indefinite Integration by Change of Variables) 


Suppose I and J are open intervals, f has domain I, g maps J into I, and g is 
differentiable on J. Assume that when we take u as the independent variable, 


[rea au = Flu) + C. 


Then when x is the independent variable and u = g(x), 
{ J(u) du = F(g(x) + C. 


PROOF Let H(x) = F(g(x)). For any x in J, the derivatives g’(x) and F’(g(x)) = (g(x) 
exist. Therefore by the Chain Rule, 


A'(x) = F'(a(x)g"(x) = S(e(x))g'(x). 
It follows that 


[ rece dx = Ho) + C = F(e(x)) 4-C. 


So when u = g(x), we have 
S(u) du = f(g(x))g'(x) dx, [re du = F(g(x)) + C. 
Theorem | gives another proof of the general power rule 


ut 1 
u" du = ——- + C, n# —l, 
| n+1 Z 
where u is given as a function of the independent variable x, from the simpler power 
rule 


xrtt 
[tax = + C¢, n#—I, 
n+l 


where x is the independent variable. 


EXAMPLE 1. Find f(4x + 1)? + (4x 4 1)? + (4x + I) dx. Let u = 4x 41. Then 
du = 4dx, dx = 1du. Hence 


ts + 18 + (4x + 1)? + (4x + Dax 


4 ue ue 


1 lfu 
= | (we +u?+u)-—d 
[ u u) qu 4 a+ 3 + 5 +C 
_ If (4x4 1)* (4x +1 (4x 4 1)? 
=i fe ge nC. 


—]1 
EXAMPLE 2 Find | —~—————~ dx. 
= l= + 1/x)* ' 


Let uw = 1 + 1/x. Then du = —1/x? dx and thus 


4.4 INTEGRATION BY CHANGE OF VARIABLES 
—1 
dx = 
x7(1 4+ 1/x)* es 


1 
du=" ~+C=-_ 1 
n= | 3 eae Celie 


= earns Pe Ter 


In a simple problem such as this example, we can save writing by using the 
term 1 + 1/x instead of introducing a new letter u, 


—1 7 1 T\.. (4 Us 
frac - Jas! +3)- —| ee 


In examples such as the above one, the trick is to find a new variable u such 
that the expression becomes simpler when we change variables. This usually must 
be done by an “educated” trial and error process. 

One must be careful to express dx in terms of du before integrating with 


respect to wu. 


EXAMPLE 3. Find f(1 + 5x)? dx. Let u= 1+ 5x. For emphasis we shall do it 
correctly and incorrectly. 


Correct: du = 5 dx, dx =}du, 
1 uw (1 + 5x)? 
1 + 5x)?d =i ee re 7 
Ja +sm x w+ du ae eee 
3 
Incorrect: fo + 5x)? dx = fe dx = 5 +C= e + +C. 
4 Z u3 (1 + 5x)? 
Incorrect: (1 + 5x)* dx = | u-du= 3 +C= 3 +C, 
EXAMPLE 4 Find { x3./2 — x? dx. Let u = 2 — x?, du = —2x dx, dx = du/(—2x). 


We try to express the integral in terms of u. 
fx aes x8 /u— . = f-5e u du. 
Since u = 2 — x? = 2 — u. Therefore 
—4x?\/u du = | —4(2 - u)/u du = { — Ju + 4? du 


= —9y3? 4 4.2y9/7 4 C 
—2(2 — x*P/2 4 42 — x”)? +.€. 


We next describe the method of definite integration by change of variables. In 
a definite integral 


i h(x) dx 


it is always understood that x is the independent variable and we are integrating 
between the limits x = a and x = b. Thus when we change to a new independent 
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variable u, we must also change the limits of integration. The theorem below will 
show that ifu = c when x = a and u = d when x = b, then c and d will be the new 
limits of integration. 


THEOREM 2 (Definite Integration by Change of Variables) 


Suppose I and J are open intervals, f is continuous and has an antiderivative 
on I, g has a continuous derivative on J, and g maps J into I. Then for any two 
points aand b in J, 

(bd) 


b 
{ CO a wee 


g(a) 


PROOF Let F be an antiderivative of f Then by Theorem |, A(x) = F(g(x)) is an 
antiderivative of h(x) = f(g(x))g’(x). Since f, g, and g’ are continuous, h is 
continuous on J. Then by the Fundamental Theorem of Calculus, 

g(b) 


b 
| S(g(x))g'(x) dx = H(b) — (a) = F(g(b)) — F(g(a)) = f(u) du. 


g(a) 


EXAMPLE 6 Find the area under the line y = 1 + 3x from x = Oto x = 1. This can 
be done either with or without a change of variables. 


Without change of variable: f(1 + 3x)dx = x + 3x?/2 + C, so 


: ax? |? 341? 3%07\. 5 
fas andv= xe 2] = (14234) [04 5 | =. 


With change of variable: Let u=1+4 3x. Then du = 3dx, dx = ¥du. 
Whenx =0,u=143-0=1.Whenx =1lu=143-1=4. 


: Ae uw|* 16 1 #15 5 
1 + 3x)dx = *—du = =—-+=—_=-, 
J (1 + 3x) dx i u 3 u | RG ETS 


Example 5 shows us that [j (1 + 3x) dx = Jf (u/3) du; that is, the areas 
shown in Figure 4.4.2 are the same. 


ypel+3x 


Figure 4.4.2 
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Figure 4.4.3 


EXAMPLE 6 Find the area under the curve y = 2x/(x? — 3)? from x = 2 to x = 3 


(Figure 4.4.3). 
Let u = x? ~— 3. Then du = 2x dx. At x =2, u=2?—3=1. At x =3, 
u = 3% — 3 = 6. Then 
7 2x 6 4 1° 1 =. 5 
bee aie Ee ed 
[ @? — 37 “* ez i Al 6° 6 


EXAMPLE 7 Find j,/1 — x? x dx. The function ,/1 — x* x as given is only defined 
on the closed interval [—1, 1]. In order to use Theorem 2, we extend it to the 


open interval J = (— 00, 00) by 
0 ifx< —-1 or x>1, 


Ho = 4 —x?x  if-1<x<1. 


Let u=1—x?. Then du = —2xdx, dx = —du/2x. Atx =0, u=1. At 


x = 1,u = 0. Therefore 


[, i= xax = [Jue bauy Si 


1 1 
= tb —1,2,,3/2 a fas 
=4/ udu = 4+4u | =3 O0=3. 
ie) 0 


We see in Figure 4.4.4 that as x increases from 0 to 1, u decreases from 1 to 0, 
so the limits become reversed. The areas shown in Figure 4.4.5 are equal. 


te) 
1 


Figure 4.4.4 
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Figure 4.4.5 
We can use integration by change of variables to derive the formula for the 


area of a circle, A = ar, where r is the radius. It is easier to work with a semicircle 
because the semicircle of radius r is just the region under the curve 


ypo,/r — x’, —rsxesr. 


To start with we need to give a rigorous definition of x. By definition, z is the area of 
a unit circle. Thus z is twice the area of the unit semicircle, which means: 


DEFINITION 
n=2f STs ax 
-1 
The area of a semicircle of radius r is the definite integral 
[ f= x? dx. 
To evaluate this integral we let x = ru. Thendx = rdu.Whenx = +r,u = +1. Thus 
[ ft? — x? dx [ Jr? = (ru)? r du = [ r?/1 — u? du 
-r -1 -1 
= ef Jl — wu? du = oS: 
-1 


Therefore the semicircle has area nr?/2 and the circle area nr? (Figure 4.4.6). 


EXAMPLE 8 Find { alee d 
—____—— qx, 
ol tafe xt 
Let w= x — x3. Then du = (1 — 3x*)dx. When x =0, u=0-O =O. 
When x = 1,u = 1 — 17 = 0. Then 


1 3x? - 1 . du 
dx | = 0. 
ol+./x—x3 o 6L + Ju 
As x goes from 0 to 1, u starts at 0, increases for a time, then drops back to 0 
(Figure 4.4.7). 
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Figure 4.4.6 


f(x) 


Figure 4.4.7 


We do not know how to find the indefinite integrals in this example. Neverthe- 
less the answer is 0 because on changing variables both limits of integration 
become the same. Using the Addition Property, we can also see that. for 
instance, 


ie 3x? —1 1 [ 3x? — 1 F 
x = X. 
2 


o 1+./x—-—x3 pil+./x—x? 


PROBLEMS FOR SECTION 4.4 


In Problems 1-90, evaluate the integral. 


1 
ee /3y +1 dy 
1 le +p? dx 2 | 3y + ldy 


3 Jo ~ 42)° dz 4 { (1 — x)? dx 


5 [2 1-Pde 6 [(—« 
f2x? +1 


4y 
22 : 
7 [xe + 5x*)’ dx 8 lo 3p dy 
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J sina ax 

{ ssincax — I)dx 
[ sine coso ao 

{ cos? @ sin 6 dO 


fe sin(x? + 1) dx 


/ sin (In x) re 


2 ¢ 


{ /sint cost dt 
i e?* dx 
[eet + be *dx 
[ xe" dx 
[be dx 


Je e*" "cos @ dO 


— Na 


ic He 4y, 
[vere 


| yJ2 + y? dy 


[-~Ra 
Jl =u? 


ds 


v/3s + 2 


|———« 
Sita ty 


[= Ge tee ha “2 dy 


10 


12 


54 


[costa — 2x)dx 

J asinx + bcosx dx 
| sin? 6 cos 6 d@ 

{ sina) + cos(30) dé 
| x? cos (x3) dx 

| e' cos(e') dt 


J s/t cos(t./t) dt 


fe l+e'dt 
2 
iE - mee 
x 
[xs i haa 
1 
lea Pim” 
1 — 2t 
fj + at 
2 
x 
Jo r 5 ax 
| 1 
——————_ dx 
(2x — 1),/1 — 2x 


See 
lars ete 


[e+ 1),/2x? + x + 5dx 


[vi — 5zdz 
1 
Ye Ea PO 
J yl — 4/y)? 


{ | 
+ dx 
/xl + 2./xy 
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59 


61 


63 


65 


67 


69 


71 


73 


75 


77 


79 


81 


83 


85 


87 


89 


feria 


[es + 6)(1 — s)-*ds 


ye 
ees 


x 
— 
es * 


f 1 — 3u du 
[= 4x -— 1 
f/4x +1 


sil a 


ye 

l+y 
6u —5 

(3u + 2)? 
J 

dx 
Pear: 


cos@ 
i sin@ a 


1 
lane 


sin@ 
1 + cos@ 


[Pea 


7 ay 


du 
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56 


38 


60 


62 


64 


66 


68 


70 


72 


74 


76 


78 


80 


82 


88 


90 


In Problems 91-108, evaluate the definite integral. 


91 


7/3 
i sin 6 d@ 


[pee 
7 
1 
if By +p?” 
2 


x 
————— a 
lee +1 * 
fev3 — x* dx 


AY 
ie + 2) as 
fv /a4 y? dy 


xs 1 
= * 
{25 Jadu 


laney + 3) o 


y 
Uu 
ire + 4) 


f A dx 

1 Eas : 

fs 4+ COSX 
———— dx 


e* + sinx 


if tan é@ dé 


(we ig 
xr+x+1 


{Ss — cos@ 
sin@ + cosé 


1 
i dx 
x Inx 


aft 
| cos(20) dé 


—nfd 
1 
f xe” dx 
0 
i 
x 
if gay rs dx 


2n 
i asin@ + bcosé dé 
0 
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2 1 
\ 
99 erin 100 { aay 
i) es 5 ax +a 


101 ic (2x + 1)? dx 102 ( t(t? + 3)~? dt 

7 2 ; 
103 [ (1 + 6x)° dx 104 i} Juri dt 

2 i x2 
105 i v./2v? + 9 dv 106 | ay® 

1 x 5 
107 i yo x 108 | x(x? + 2)'3 dx 

‘ ; Jé 

109 Find the area of the region below the curve y = 1/(10 — 3x) from x = 1 to x =2. 
110 Find the area of the region under one arch of the curve y = sin x cos x. 
111 Find the area of the region under one arch of the curve y = cos (3x). 
112 Find the area of the region below the curve y = 4x./4 — x? between x = Oand x = 2. 
113 Find the area below the curve y = (1 + 7x)?/3 between x = Oand x = 1. 
114 Find the area below the curve y = x/(x? + 1)? between x = O and x = 3. 


115 Evaluate: [ eee Nac ee A oily 
“10 1 ee 


1 
116 Evaluate: i} 2x./(1 — x?)> + 1 dx 
-1 


117 Let f and g have continuous derivatives and evaluate | f“(g(x))g'(x) dx. 

118 A real function f is said to be even if f(x) = f(—x) for all x. Show that if f is a continuous 
even function, then f°, f(x) dx = J? f(x) dx. 

119 An odd function is a real function g such that g(—x) = — g(x) for all x. Prove that for 


a continuous odd function g, f°, g(x) dx = 0. 


AREA BETWEEN TWO CURVES 


A region in the plane can often be represented as the region between two curves. 
For example, the unit circle is the region between the curves 


y=—f/l—x*, y= /l—x?, -lexsl 


shown in Figure 4.5.1. Consider two continuous functions fand g on [a, b] such that 
I(x) < g(x) for all x in [a,b]. The region R, bounded by the curves 


y=f@), y=g%), *x*x=a x=6, 


is called the region between f(x) and g(x) from a to b. If both curves are above the 
x-axis as in Figure 4.5.2, the area of the region R can be found by subtracting the 
area below ffrom the area below g: 


b b 
area of R = i g(x) dx — | T(x) dx. 
It is usually easier to work with a single integral and write 


b 
area of R = I (g(x) — f(x)) dx. 
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y 


f) 


| 
I 
| 
! 
| 
1 
I 
| 
| 
| 
| 
| 
| 
i 


“y Seaeee 
of 


Figure 4.5.1 Figure 4.5.2 


In the general case shown in Figure 4.5.3, we may move the region R above the 
x-axis by adding a constant c to both f(x) and g(x) without changing the area, and 
the same formula holds: 
b b 
area of R = | (g(x) + e) dx — Kero) + dx 


b 
= i (a(x) — fx) dx. 


y 


y@) +e 


fx) 


Figure 4.5.3 


To sum up, we define the area between two curves as follows. 
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DEFINITION 


If f and g are continuous and {(x) < g(x) for a < x S b, then the area of the 
region R between {(x) and g(x) from a to b is defined as 


b 
I (g(x) — f(x)) dx. 


EXAMPLE 1. Find the area of the region between the curves y = 4x* — l and y = x 
from x = 1 to x = 2. In Figure 4.5.4, we sketch the curves to check that 
4x? — 1 <x for! <x <2. Then 


2 
A= ie (4x? Dax = ft be +3] = 8, 
i 


1 


Figure 4.5.4 


EXAMPLE 2. Find the area of the region bounded above by y = x + 2 and below 
by y = x?. 
Part of the problem is to find the limits of integration. First draw a sketch 
(Figure 4.5.5). The curves intersect at two points, which can be found by 
solving the equation x + 2 = x? for x. 


x? — (x + 2)=0, (x + 1)(x — 2) =0, 


x=-I1 and x=2. 
2 


2 
Then A=| (+2 — deed? +2 he] ee 
=e 


~1 


Figure 4.5.5 
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EXAMPLE 3 Find the area of the region R bounded below by the line y = —1 and 
above by the curves y = x? and y = 2 — x. The region is shown in Figure 
4.5.6. 


Figure 4.5.6 


This problem can be solved in three ways. Each solution illustrates a different 
trick which is useful in other area problems. The three corners of the region 


are: 
(—1, —1), where y = x? and y = —1 cross. 
(3, — 1), where y = 2 — x and y = —I cross. 
(1, 1), where y = x? and y = 2 — x cross. 


Note that y = x? and y = 2 — x can cross at only one point because x? is 
always increasing and 2 — x is always decreasing. 


FIRST SOLUTION Break the region into the two parts shown in Figure 4.5.7: 
R, from x = —1 tox = 1, and R, from x = 1 tox = 3. Then 


area of R = area of R, + area of Rp. 
1 1 
area of R, = | 8 — (Dax = st + | im 
-1 —] 
3 


3 
area of R = [ (2 — x) —( Idx = 3x — br] = 2. 
1 


1 


areaofR=242=4. 


Figure 4.5.7 First solution 
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SECOND SOLUTION Form the triangular region S between y = —1 and y=2—x 
from —1 to 3. The region R is obtained by subtracting from S the region S, 
shown in Figure 4.5.8. Then 


area of R = area of S — area of S,. 


3 3 
area of S = | (2 — x) — (—1)dx = 3x i = 8. 
=I 


1 if 
area of S, = | (2 — x) — x3. dx = 2x — 4x? is = 4, 
raul | 


areaofR=8 —4=4, 


Figure 4.5.8 Second solution 


THIRD SOLUTION Use y as the independent variable and x as the dependent 
variable. Write the boundary curves with x as a function of y, 


yp=2- x becomes x = 2 — y. 
ee becomes x = y!%, 
The limits of integration are y = —1 and y = 1 (see Figure 4.5.9). Then 


1 1 
A= | (2 — y) — yx"? dy = 2y — dy? — ayn] = 4, 
-1 


~1 


As expected, all three solutions gave the same answer. 


¥ 


{ \ 


aaa | 
Figure 4.5.9 | Third solution 
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PROBLEMS FOR SECTION 4.5 


In Problems 1-43 below, sketch the given curves and find the area of the region bounded by 


them. 


Oo Pert Aun bk WN = 


Oe oe a 
SSCarAA AWN KK S 


Ne 
—_ 


r 8 BR 


® yy NN bY NV 
— | 


f(x) =0, g(x) =5x—x?, OXx<4 
fe) =/x%, x) =x4, 1S x54 


fx) =x/1—x?, gx) =1, -1 Sxl 


yp=x—-2, y=3x'?, O<x<l 
y=4/x, y= /xt+1, OSx<4 


y= fe+l—x yoJ/x*?t1t+x —-0Sxs1 


The x-axis and the curve y = —5 + 6x — x 


The x-axis and the curve y = 1 — x* 

The y-axis and the curve x = 25 — y* 

The y-axis and the curve x = y(8 — y) 
y=cosx, y=2cosx, —2z/2<x< 2/2 
y=sinxcosx, y=1l, O<Sx<n/2 
y=-sinx, yp=sinx, 0OSx<n 
y=sinx, y=cosx, 0Sx< 2/4 


y=sinxcosx, p=sinx, O<x<n 


y=sin?xcosx, y=sinxcosx, 0<x<2/2 


y=x, py=e, OSxS2 
y=e*, p=e, OSx<2 
y= -e%, y=e, -1sx<i 
y=xe", p=e O0<x<l 


1 


y y=1, O<Sx<2 


“547° 
1 1 
? De es eee Seo 
y=1/x, y=x, 1Sx<2 
y=, y=h, O<xS1 
x* +1 


f(x) = x97, g(x) = x? 
yp=x?—2x, p=x—2 
ypaxt—2x?, p= 2x? + 12 
y=xt-1, p=xP—-x 
y=x4o2 +0, y=1/e? +) 


y=x? /1 x, y=x 1 — x’, O<x 


y=2x?, y=x?4+4 


x=y, x=2-y 

Jx + Jy = 1 and the x- and y-axes 
x?y=4, x? +y=5 (first quadrant) 
yp=x/xt1, y= 2x 

y=0, y=xetxt2, x=2 
y=2x+4, y=2—3x, y= -x 
y=x?-1, y=(x- 1), y=(x4+ 1? 
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39 y Js y=, p= 10- 2x 

40 ySex=h pedo, yea/x 

41 ye-x, ypoJx, y=3x-2 

42 y= -2, p=ortx xty=3 

43 yor’, p=2x-2, y= 2x73 (first quadrant) 

44 Find the area of the ellipse x?/a? + y?/b? = 1. Use the fact that the unit circle has area x. 

45 Sketch the four-sided region bounded by the lines y= 1, p=x, = 2x, and 
y = 6 — x and find its area. 

46 Find the number c > O such that the region bounded by the curves y = x, = —2x,and 
x = c has area 6. 

47 Find the number c > 1 such that the region bounded by the curves y = 1, y = x77, 
and x = c has area 1. 

48 Find the number ¢ such that the region bounded by the curves y = x? and y' = c has 
area 36. 

49 Find the number c > 0 such that the region bounded by the curves y = x* and y = ex 
has area 9. 

50 Find the value of ¢ between —1 and 2 such that the area of the region bounded by the 
lines y = —x, y = 2x, and y = 1 + cx isa minimum. 

51 Find the value of c such that the line y = c bisects the region bounded by the curves 


y=x?andy=1. 
§2 Find the value of ¢ such that the line y = cx bisects the region bounded by the x-axis 
and the curve y = x — x?. 


NUMERICAL INTEGRATION 


In numerical integration, one computes an approximate value for the definite 
integral rather than finding an exact value. In this section we shall present two 
methods of numerical integration, called the Trapezoidal Rule and Simpson’s Rule. 

The Fundamental Theorem of Calculus gives us a method of computing 
the definite integral of a given continuous function f from a to b. The method is to 
find, by trial and error, an antiderivative F of fand then to use the equation 


b 
i) f(t) dt = F(b) — F(a). 


When the method works, it provides an exact value for the integral. However, the 
method succeeds only if the antiderivative happens to be a function that can be 
described in a simple way. For many integrals one cannot find a formula for the 
antiderivative, and the method fails. Such integrals can still be computed approxi- 
mately using numerical integration. 

The Trapezoidal Rule and Simpson’s Rule can always be applied and do 
not use the antiderivative. They are easy to carry out on acomputer or hand calculator. 
We already discussed one method of approximating the definite integral in Section 4.1, 
the Riemann sum. The Trapezoidal Rule is a modified form of the Riemann sum, 
which gives a much closer approximation for a given amount of effort. Simpson’s 
Rule is a further modification that gives still better approximations. 
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Let f be a continuous function on an interval J, and let a < b in I. By 
definition, for each positive infinitesimal dx the definite integral 


ie f(x) dx 


is the standard part of the infinite Riemann sum 
b 
d f(x) dx, 
b b 
{ f(x) dx = ofS f@) ax. 


In Section 4.1, examples were worked out to show that the finite Riemann sums 
become very close to the definite integral when Ax is small; that is, the finite Riemann 
sums approximate the definite integral. In Section 4.2, we saw that the definite 
integral is the limit of the finite Riemann sums as Ax > 0°: 


b b 
i f(x) dx = lim 3 f(x) Ax. 


The Riemann sum, which is a sum of areas of rectangles, is a rather inefficient 
approximation of the definite integral. We can usually get a much closer approxi- 
mation with the same amount of work by adding up areas of trapezoids instead of 
rectangles, forming the Trapezoidal Rule suggested by Figure 4.6.1. The Trapezoidal 
Rule also provides a formula, called an error estimate, which tells us how close the 
approximation is to the exact value of the definite integral. 


Riemann Sum Trapezoidal Approximation 


Figure 4.6.1 


Choose a positive integer n and divide the interval (a, b] into n subintervals 
of equal length Ax = (6 — a)/n. The partition points are a = X9,X1,X2,-..,X, = b. 
The trapezoidal approximation is the area of the region under the broken line con- 
necting the points 
(xo, f (Xp), (x1 » F(x), eee (Xn>f (Xq))- 
Since all of these points lie on the curve y = f(x), the broken line closely follows the 


curve. So one would expect the area of the region under the broken line to closely 


approximate the area under the curve. 
Consider a single subinterval [x,,, X,,4,] of width Ax. The region under the 
line segment connecting the two points 


ings F (Xm) (Xin+ 1 st ne 1)) 
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is a trapezoid and its area is 


Sn) +O ns 1) AD 


The sum of the areas of the trapezoids is a modified Riemann sum 


> f(x) + f(x + Ax) 
> 3 A 


a 2 


7 I (Xo) +1) 4 (x4) + Le) Beicett et ee Re 


= [Ff (Xo) + f(%1) + fla) Fo + S$ %n=1) + 3S (%,)] AX. 
We thus make the definition: 


DEFINITION 


Let Ax = (b —a)/n evenly divide b — a. Then by the trapezoidal approxi- 
mation to the definite integral [? f (x) dx we mean the sum 


pete 


Ax = 5 (Xo) a J (x1) a ae f(%,-1) + sf (x,)] Ax. 


The Trapezoidal Approximation of an integral ae ge x)dx can be computed 
very efficiently on most hand calculators. First compute the sum 


3f (Xo) + f(X1) + f(x.) ++ + 4F(%,) 


by cumulative addition. Then multiply this sum by Ax to obtain the Trapezoidal 
Approximation. 


THEOREM 1 


For a continuous function f on (a, b], the trapezoidal approximation approaches 
the definite integral as Ax > 0*, that is, 


5 £9) + f(t AX) Ay 


b 
i {[(x)dx = lim 
a Ax> Ot 2 
PROOF Comparing the formulas for the trapezoidal approximation and the Riemann 
sum, we see that 
> f(x) + f(x + Ax 2 
yO FEST ES ax =F fla) Ax + ES) — EF 60)) As 


a 


For dx positive infinitesimal, the extra term 
(of %u) — f(x) dx 
is infinitely small. It follows that 


b Lf, 9 : 
SLO TIE TO) wey sexydy = [foods 
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From a practical standpoint, it is desirable to have a good estimate of error. 


We shall first work an example and then state a theorem which gives an error estimate 
for the trapezoidal approximation. 


EXAMPLE 1 Approximate the definite integral 


1 
{ /1 + x? dx. 
oO 


Use the trapezoidal approximation with Ax = 4. We first make a table of 
values of ./1 + x*. The graph is drawn in Figure 4.6.2. 


JESU 8 dx ~ 1.1501 


x 
Figure 4.6.2 
¥ ee | Jl+x?to | term in trapezoidal 
four places approximation 
Xp =0 1 1.0000 0.5000 = $ f(x) 
x, =4 ./1.04 1.0198 1.0198 = f(x,) 
x, =2 1.16 1.0770 1.0770 = f(x) 
x,=3 1.36 1.1662 1.1662 = f(x) 
<4 1.64 1.2806 1.2806 = f (x4) 
x, =1 <f2 1.4142 0.7071 = 4 f (xs) 
5.7507 = total 


Thus, 3f (xo) + f(x1) +f.) + f(%3) + fa) + 3/5) = 5.7507. Since 
Ax = 4, the trapezoidal approximation is 


(5.7507) * 4 = 1.1501, 


1 
{ Jl +x? dx ~ 1.1501. 
0 


The trapezoidal approximation can be made as close to the definite integral 
as we want by taking Ax small. From a practical standpoint, however, it is helpful 
to know how small we should take Ax in order to be sure of a given degree of accuracy. 
For instance, suppose we need to know the definite integral to three decimal places. 
How small must we take Ax in our trapezoidal approximation? The answer is given 
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by the Trapezoidal Rule, which gives an error estimate for the trapezoidal 
approximation. 
The error in the trapezoidal approximation is the absolute value of the 
difference between the trapezoidal sum and the definite integral, 
b x A b 
error =|) £@) +f =hy) Ax i} f(x) dx 


a 


An error estimate for the trapezoidal approximation is a function E(Ax), which is 
known to be greater than or equal to the error. 

Thus if E(Ax) is an error estimate, the trapezoidal sum is within E(Ax) of 
the definite integral. If we want to be sure that the trapezoidal approximation is 
accurate to three decimal places—i.e., the error is less than 0.0005—we choose 
Ax so that E(Ax) < 0.0005. We are now ready to state the Trapezoidal Rule. 


TRAPEZOIDAL RULE 
Let f be a function whose second derivative f" exists and has absolute value 
at most M on a closed interval [a, b], 
\f’C)| <M fora<x <b. 


If Ax evenly divides b — a, then the trapezoidal approximation of the definite 
integral of f has the error estimate 


b-a 5 


That is, 


b fry : . b 
5 LOFTE FAY yf peayde| <P! Mary. 


a 


The proof is omitted. 


EXAMPLE 1 (Concluded) Welet f(x) = ./1 + x*. Then 


; x 
(x = I? 
Tice 
("x _ ft txt = x7/J/1 t+ x? | 
i x) I+ x? = (1 =p. x?)3/2' 


Therefore | f"(x)| < 1 for all x in [0,1]. We take M = | and use the error 
estimate given by the Trapezoidal Rule, 


b-—a if :\7 1 
M z= cate a =. . 
iG = (;| 3 


Thus our approximation is within an accuracy of 1/300, 


bg he | 
{ J1 + x? dx — 1.159 < 1/300 ~ 0.0033. 
is) 


This shows that the integral is, at least, between 1.146 and 1.154. 
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In this particular example we can even conclude that the integral is between 
1.146 and 1.150 (rounded off to three places). That is, the integral is less than its 
trapezoidal approximation. This is because the second derivative f"(x) = (1 + x*)~ 3? 
is always greater than 0, whence the curve is concave upwards and therefore y = f(x) 
is always less than or equal to the broken line used in the trapezoidal approximation. 
Actually, the value to three places is 1.148. This can be found by taking Ax = 75. 


EXAMPLE 2 Consider the integral 


i ./1 — x? dx = 2/2. 
Let fi) = /1 — x?. 


By Theorem 1, we have 


1 
; Six) + f(x + Ax), 
Paci 2 BN 5F 


However, the Trapezoidal Rule fails to give an error estimate in this case 
because f’(x) is discontinuous at x = +1. 


We now turn to Simpson’s Rule, for which the number of subintervals n 
must be even. As before, we divide the interval [a,b] into subintervals of equal 
length Ax with the n + 1 partition points 


a= Xo, Xq5--+5Xy = b, 
We shall use subintervals of length 2 Ax rather than Ax. On each of the n/2 sub- 
intervals 
[Xo; X2], [x2, Xa], Laces ) [x,=25 Xuls 


of length 2 Ax we approximate the curve y = f(x) by a parabolic arc that meets the 
curve at both endpoints and the midpoint of the subinterval, as shown in Figure 4.6.3. 
We then add up the areas under each of the parabolic arcs to obtain an approximation 
to the area under the curve, which is the definite integral. We begin with a lemma that 
gives a formula for the area of the region under one parabolic are. 


Figure 4.6.3 


LEMMA 


The area of the region under the parabola through three points (u, r), (u + h,s), 
and (u + 2h, t) (shown in Figure 4.6.4) is 
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230 4 INTEGRATION 


(ur) 


(a +2h,0 


(a +h,s) 


Figure 4.6.4 u uth u + 2h 


} 
xG ace), 


The lemma is proved at the end of this section. Using the lemma, we find 
that the area of the region under one parabolic arc from x, to x, 4.2 is 


Ax 
> Loe) + 46 ear) +f Oued]. 


It follows that the sum of the n/2 regions under the parabolic arcs is a modified 
Riemann sum, 


3 


aM|> 


- + 4f(x + Ax) +f(x + 2 =) as 


A 
= = {Lf (xo) + 4f (1) + fa] + Lf G62) + 4f 3) +f e4)] £23 


Ax 
= J [1 Co) + 41) + Fa) + Os) + Ora) to + FC) + FO) 
This modified Riemann sum is Simpson’s approximation to the definite 
integral. Note the sequence of coefficients, 
1,4,2,4,2,...,2,4, 1. 


Like the trapezoidal approximation, it is easily computed on a computer or hand 
calculator. 


THEOREM 2 


For a continuous function f on [a,b], Simpson’s approximation approaches 
the definite integral as Ax > 07, 


? A 
[foo dx = tim, S £6 €%0) + Afr) + 2F 6) + 4F 3) + FC6N 


Simpson’s approximation is almost as easy to calculate as the trapezoidal 
approximation, but is much more accurate. Simpson’s Rule is an error estimate 
that involves the fourth derivative of the function and the fourth power of Ax. 
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SIMPSON‘S RULE 
Suppose the function f has a fourth derivative on the interval [a,b] that has 
absolute value at most M, 
|fR)| < Mfora<x<b. 


If (a, b] is divided into an even number of subintervals of length Ax, then 
Simpson’s approximation to the definite integral has the error estimate 
b—~a 
180 


M(Ax)*. 


EXAMPLE 3 Use Simpson’s Rule with Ax = 0.25 to approximate the integral 
1 
A= [ e7* 2 dy 
0 


and find the error estimate. 

The curve is the normal (bell-shaped) curve used in statistics, shown in 
Figure 4.6.5. We are to divide the interval [0,1] into four subintervals of 
equal length Ax = 0.25. The following table shows the values of x and y 
and the coefficient to be used in Simpson’s approximation for each partition 
point. 


Figure 4.6.5 Example 3 


x en l2 Coefficient 
0.0 1.000000 1 
0.25 0.969233 4 
0.5 0.882496 2 
0.75 0.754840 4 
1.0 0.606531 1 


The sum used in the Simpson approximation is then 


[1.000000 + 4 + (0.969233) + 2 « (0.882496) + 4+ (0.754840) + 0.606531] 
= 10.267816 


PROOF 


4 INTEGRATION 


To get the Simpson approximation, we multiply this sum by Ax/3: 
S = (10.267816) + (0.25)/3 = 0.855651. 


To find the error estimate we need the fourth derivative of 


-x7/2 


ye 
The fourth derivative can be computed as usual and turns out to be 
yh = (x* = 6x? ae Sie <2. 


On the interval [0,1], } is decreasing because both x* — 6x? + 3 and 
—x?/2 are decreasing, and therefore y‘* has its maximum value at x = 0 
and its minimum value at x = 1, 


maximum: —_—y“#(0) = 3 
minimum: y1) = — 1.213061 


The maximum value of the absolute value | y“| is thus M = 3. The error 
estimate in Simpson’s Rule is then 
b-—a 1 
2S (Ax)*M = 
iso “4*) 180 


This shows that the integral is within 0.000065 of the approximation; that is, 


+ (0.25)* «3 = 0.000065. 


1 
[ e 7/2 dx = 0.855651 + 0.000065, 
) 


or using inequalities, 
1 


0.855586 < { el? dx < 0.855716. 


ie) 


For comparison, a more accurate computation with a smaller Ax shows 
that the actual value to six places is 


1 
i el? dx = 0.855624. 
) 


The Trapezoidal Rule for this integral and the same value of Ax = 0.25 
give an approximate value of 0.85246 for the integral and an error estimate 
of 0.00521. 


OF THE LEMMA _ The algebra is simpler if the y-axis is drawn through the 
second point, so that u +h = 0, and the three points have coordinates 


(—A, r), (0, s), (A, 0). 


Suppose the parabola has the equation y = ax? + bx +c. Then the area 
under the parabola is 


h 
A= | (ax? + bx + c)dx 
—h 


oe DS aires 
| 2 


-h 


2 3 
=e ah? + 2ch. 
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When we substitute the coordinates of the three points (—4,r), (0, 5), (h, 0) 
into the equation for the parabola, we obtain the three equations 


r = ah? ~ bh+c, 
S=C, 
t=ah? + bh+e. 
Add the first and third equations and solve for a: 
r +t = 2ah? + 2c 
ies tea 
“eh 


Finally, substitute the above expression for a and s for c in the equation 
for the area: 


2 
A= zane + 2ch 


r+t~—2e 2 ,, 
= yer eae + 2ch 
r+t—2¢+ 6c 
Ree eae = oe 


=F + 4c + 1). 


h 
=a +4st 0). 


PROBLEMS FOR SECTION 4.6 


Approximate the integrals in Problems 1~20 using (a) the Trapezoidal Rule and (b) Simpson’s 
Rule. When possible, find error estimates. If a hand calculator is available, do the problems 
again with Ax = 0.1. 


1 


11 


13 


2 2 
{ xax, Ax = 0.5 2 { x°ax Ax = 0.5 
0 0 
2 3 4 
( x? —1dx, Ax = 0.25 4 f <a Ax = 0.5 
‘ 1x 
2 L 1 
| zdx, Ax = 0.25 6 [x x+1dx, Ax = 0.25 
i i+x 0 
5 1 
{ je Ax = 0.5 8 [ Ve ¥ Tax, Ax =4 
1x41 o 
1 4 
i xt 4+1dx, Ax =} 10 i 14+ 1/xdx, Ax =0.5 
0 1 
6 | 12 j 
| dx, Ax=1 12 i) dx, Ax =2 
o xt] 9 2x+3 
13 1 4 1 
i Sa gy ies 14 { dx, Ax =1 
i eee @ 24k/k 


233 


234 


15 


17 


19 


21 


22 


23 


24 


25 


EXTRA 
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i sin@dé, Ax =n/2, 2/10 16 | sin? @d0, Ax =n/2, n/10 
0 0 
1 1 
{ e“dx, Ax=+ 18 i) e dx, Ax =4 
0 0 
2 2 
I Inxdx, Ax =4 20 { In(i/x)dx, Ax =4+ 
1 1 


Let f be continuous on the interval [a, b] and let Ax = (b — a)/n where n is a positive 
integer. Prove that the trapezoidal sum is equal to the Riemann sum plus 4{ f(b) — 
I (a)) Ax, that is, 

b 


b 
LEU) + f(x + Ax) Ax -[d 09 Ax 


a 


+ 3( f(b) — f(a) Ax. 


Show that if f(a) = f(b) then the trapezoidal sum and Riemann sum are equal. 


Prove that for a linear function f(x) = kx + c, the trapezoidal sum is exactly equal to 
the integral. 


Show that if f(x) is concave downward, f(x) > 0, then the trapezoidal sum is less 
than the definite integral of f(x). 


Show that for a quadratic function f(x) = ax? + bx + c, Simpson’s approximation 
is equal to the definite integral. 


Show that for a cubic function f(x) = ax? + bx? + ex + d, Simpson’s approximation 
is still equal to the definite integral. 


PROBLEMS FOR CHAPTER 4 
Evaluate Tid, Ax = 1/4 


Evaluate Ties Ax, Ax =2 


Evaluate )? ,2*Ax, Ax =1 


Evaluate 35 x,/x + 1 Ax, Ax = 1/2 

If F(x) = 1/(2x — 1)? for all x # 1/2, find F(2) — F(1). 

If G(t) = ./4t + 1 for allt > —1/4, find G(2) — G(0). 

A particle moves with velocity v = (3 + 2/1). How far does it move from times ty = 1 
tot, = 5? 

A particle moves with velocity v = t?,/t? — 1. How far does it move from times ty = 1 
tot, = 4? 

A particle moves with velocity v = (f + 1)(2r + 3). If it has position y) = 0 at time 
t = 0, find its position at time t = 10. 


A particle moves with acceleration a = 1/t*. Ifit has velocity »p = 4and position yy = 2 
at time t = |, find its position at time t = 3. 


Find the area of the region under the curve y = oe. l<x<4 
Find the area of the region under the curve y = Jx — x./x, O<x<l. 


In Problems 13-30, evaluate the integral. 


13 


[a -2@ + 3yax 14 Je+3e-2 re 


Oo 


40 


41 


42 


43 


44 
45 
46 


47 


EXTRA PROBLEMS FOR CHAPTER 4 


| aoe 16 | (4x + 1)! dx 

| (u/,/1 — 3u?) du 18 (agree dx 
[oat Ha / de Dat 20 a ee 
[ovr¥2ey 2 Ja — fx)~+ dx 
{ cos (5) dx 24 [ vesin Jx dx 

{ etd 26 | te F 


4 6 
{ (y + /y)dy 28 i (x/./x? — 1) dx 
0 2 


1 1 
| e** dx 30 | x sin (x?) dx 
0 


iv) 


0 


4 y 
Differentiate [ wh Be /x — 1dx 34 Differentiate [ (If(x + ./x) dx 


y 


x 3x 
Differentiate [ Jt + 2dt 32 Differentiate | (t2/(t2 — 1))dt 
1 


Find the function F such that F(x) = x — 1 for all x, and the minimum value of F(x) is b.. 
Find the function F such that F’(x) = x for all x, F(0) = 1, and F(1) = 1. 
Find the function F such that F’(x) = 6 for all x, F(x) has a minimum at x = 1, and the 
minimum value is 2. 
Find all functions F such that F’(x) = 1 + x7? for all positive x. 
Find the function F such that 
‘ 0 ifx <0 

POs ” ifx >0 
and F(0) = 1. 
Find the value of b such that the area of the region under the curve y = x(b — x), 
0<x <5b,is1. 
Suppose fis increasing for a < x < b, and Ax = (b — a)/n where nis a positive integer. 
Show that 


< [f@) — f(@] Ax 


5 b 
Y f(x) Ax — | roves 
Suppose f is continuous for a < x < b. Show that 
b b 
[ senax] < [ rcatax. 
Find the area of the top half of the ellipse x*/a? + y?/b? = 1 using the formula 
n= 21, /1 — uw? du. 
Evaluate f* , (1 — x)*(1 + x)’ dx using the formula x = 2f1, ./1 — u? du. 
Find dy/dx if y = Ss x f(t) dt. 
Suppose f(t) is continuous for all ¢ and let G(x) = f(x — t)f(t)dt. Prove that 


G"(x) = f(%). 


Prove that for any continuous functions f and g, 
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b 


b 
Sf (x)e(x) dx < af f(x dx + Bf a(x)? dx. 


ab 
2AB | 


a 


48 Prove Schwartz’ Inequality, 


b b- “ab 
[ Txdglx) dx < | { I(x as | g(x)? dx. 


va 


Hint: Use the preceding problem. 
OO 49 Suppose f is continuous and dx is positive infinitesimal. Show that 


b b 
Y f(x + Zdxjdx = f(x) dx. 
Hint: For each positive real c, 
1 
f(x) -—e¢ <1(x + 54x] <f(x) te. 


Use this to show that 


b b 1 b 
[ {[Q)dx —cb-a<¥ f (x + 54x) dx < | I(x) dx + cb — a). 


’ 


O 50 Suppose f is continuous, n is an integer, and dx is positive infinitesimal. Prove that 


b b 
Y fle + ndx) dx = { f(x) dx. 


